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Abstract - An analysis is made for transient 
Magnetohydrodynamic (MHD) flow of two immiscible fluids 
thorough a horizontal channel. The flow above the interface is 
assumed to be electrically conducting while the other fluid and 
the channel walls are assumed to be electrically insulating. The 
partial differential equations governing the flow are solved 
analytically using two-term periodic and non-periodic functions 
in both the regions of the channels. Velocity profiles are 
analyzed for various physical parameters such as viscosity ratio, 
frequency parameter and Hartmann number. 
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I. INTRODUCTION 
Unsteady flow of immiscible fluids through a channel is 
of special importance in the petroleum extraction and transport, 
the reservoir rock of an oil field always contains several 
immiscible fluids in its pores. Lighthil [1] was the first to have 
studied the unsteady forced flow of a viscous incompressible fluid 
past a flat plate and a circular cylinder with small amplitude 
oscillation in free stream. The corresponding problem of unsteady 
free convection flow along a vertical plate with oscillating surface 
temperature was studied by Nanda and Sharma [2]. Later Muhuri 
and Maiti [3] and Verma [4] have analyzed the effect of 
oscillation of the surface temperature on the unsteady free 
convection from a horizontal pipe. Hossain et al [5] used the 
linearised theory to study the free convection boundary layer flow 
of electrically conducting fluid along a vertical plate to surface 
temperature oscillations. 

All these studies pertain to one fluid model. Most of the 
problems relating to petroleum industry, geophysics, plasma 
physics, magneto-fluid dynamics etc, involve multi-fluid flow 
model. Both theoretical and experimental work is found in 
literature on stratified laminar flow of two immiscible fluids in a 
horizontal pipe. In the petroleum industry, as well as in other 
engineering fields, stratified two-phase flow often occurs. There 
have been some experimental and analytical studies on 
hydrodynamic aspects of two-phase flow reported in the recent 
literature. The interest in these types of problems stems from the 
possibility of reducing the power required to pump oil in a 
pipeline by suitable addition of water. The first investigations 
were associated with the LM-MFM generator project at the 
Argone National Laboratory. Rudraiah et al [6] studied the 
Hartmann flow past a permeable bed in the presence of a 
transverse magnetic field with an interface at the surface of the 
permeable bed. Packham and Shail [7] analyzed stratified laminar 
flow of two immiscible liquids in a horizontal pipe. Hartmann 
flow of conducting fluid in a channel with a layer of non- 



conducting fluid between upper channel wall and the 
conducting fluid has been studied by Shail [8]. He found 
that increase of order 30% could be achieved in the flow 
rate for suitable ratios of depths and viscosities of two 
fluids. Loharsbi and Sahai [9] dealt with two-phase MHD 
flow and heat transfer in a parallel-plate channel. Both 
phases are incompressible and the flow was assumed to be 
steady, one-dimensional and fully developed. The study 
was expected to be useful in the understanding of the 
effect of the presence of slag layers on the heat transfer 
characteristics of a coal-fired MHD generator. Alireza and 
Sahai [10] studied the effect of temperature -dependent 
transport properties on the developing MHD flow and heat 
transfer in a parallel plate channel whose walls are held at 
constant and equal temperatures. Following the work of 
Alireaz and Sahai [10], Malashetty and Umavathi [11] and 
Malashetty et al [12-14] studied two-phase MHD flow and 
heat transfer in an inclined channel. Stamenkovic,M. Z., 
et. al [15] investigated the MHD flow of two immiscible, 
electrically conducting fluids between isothermal and 
insulated moving plates in the pressure of an applied 
electric and inclined magnetic field with the effect of 
induced magnetic field. 

The above investigations are related to that the 
flow is steady. However the study closed to practical 
problems is to consider unsteadiness. Chamkha [16] 
studied unsteady MHD convective heat and mass transfer 
past a semi-infinite vertical permeable moving plate with 
heat absorption. Recently Umavathi et al [17-19] have 
presented analytical solutions for unsteady/oscillatory two- 
fluid flow and heat transfer in a horizontal channel. 
Keeping in view the wide range of applications in 
geophysics and MHD generators, an attempt is made to 
analyze flow nature and heat transfer of unsteady 
Hartmann flow of two immiscible fluids in a horizontal 
channel. 

II.MATHEMATICAL FORMULATION 

The geometry under consideration consists of 
two infinite parallel plates extending in the Z and X 
directions. The regions < y < h and -h < y < are 
denoted as Region-I and Region-II respectively. The 
fluid in Region-I is a conducting fluid having density p, 

and viscosity ju l . A constant magnetic field of strength B 
is applied transverse to the flow field. The Region-II is 
filled with a non-conducting fluid having density p 2 and 

viscosity ju 2 ■ It is assumed that the flow is unsteady, fully 
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developed and that all fluid properties are constants. The flow is 



assumed to be driven by a pressure gradient 



V 



dx 



for region-I 



and 



Sp 2 
dx 



for region-II. 



u l (h)=0 
u 2 (- h) = 

M 1 (0) = M 2 (0) 



9m, 



a 



du 



dy dy 



(7) 



aty = 



Under these assumptions and taking p x = p 2 = p the 

governing equations of motion and energy are given by: 
Region-I 

du l 



Po 



dt 



■7 L + ATy-Wi 

ox dy 



Region-II 



Po 



du 2 
~d7 



dx dy 2 



(1) 



(2) 



where u is the x-component of fluid velocity. The subscripts 1 and 
2 correspond to region-I and region-II, respectively. The 
boundary conditions on velocity are the no-slip boundary 
conditions which require that the x-component of velocity must 
vanish at the wall. We also assume continuity of velocity and 
shear stress at the interface between the two fluid layers at y=0. 

The hydrodynamic boundary and interface conditions 
for the two fluids can then be written as 
u^h) = 
u 2 (-h)=0 

u l (o) = u 2 (o) (3) 

du, du ? 
Mi — = Mi— aty = 

dy dy 



III. SOLUTIONS 
The governing momentum equations (5) and (6) 
are solved subject to the boundary and interface conditions 
(7) for the velocity distributions in both regions. These 
equations are partial differential equations that can not be 
solved in closed form. However, it can be reduced to 
ordinary differential equations by assuming 

u i (y,t)=u i0 +e Uot u a{y) (8) 



P ! .(y,f)=/ , ! o+e'' ffl ' J P, 1 ( > ) 
fori = 1,2 



(9) 



By substitution of equation (8) and (9) into 
equations (5) and (6), one obtains the following pairs of 
equations: 

Region-I 

Non-periodic coefficients 



d 2 u 



10 



dy 1 



-M L u 10 +P w =0 



Periodic coefficients 



d m, 



dy 



^--{M 2 + ico) u n + P n =0 



(10) 



(11) 



By using the following non-dimensional quantities: 

* h 2 * 

u i =U a u i , y = hy , t = — t , 
v 



hi 



dx 



M = B Q h 



(4) 



/Uy U Q \ dX J \ //[ 

and for simplicity dropping the asterisks, equations (1) and (2) 

become 

Region-I 



du, n d 2 u, ,0 
dt 1 dy 2 



Region-II 



du 2 
~dt 



P 2 +a 



d 2 u 7 

~d7 



u 2 

a = — is the ratio of viscosity. 
Mi 



(5) 



(6) 



Region -II 

Non-periodic coefficients 



a ■ 



,2 

d u 



dy 



2 - + ^20=0 



Periodic coefficients 



d 2 Un 



a- 



dy 2 



■ ico u 2l + P 2l =0 



(12) 



(13) 



The corresponding boundary and interface 
conditions become as follows: 
Non-periodic coefficients 
u l0 {h)=0 

u 2Q (-h)=0 

«io(o)=«2o(o) (14) 



du 



dy 



10 _ 



= a- 



du 



2(1 



dy 



aty = 



The hydrodynamic boundary and interface conditions for both 
fluids in non-dimensional form become 



Periodic coefficients 
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W n (/z) = 

u 2 i(-h) = 
«n(0) = «2i(0) 



(15) 



a- 



du 



21 



dy 



aty=0 



Equations (10) to (13) along with boundary and interface 
conditions (14) and (15) represents a system of ordinary 
differential equations and conditions that can be solved in closed 
form. Since the solutions can be obtained directly and depicted 
graphically. 

HI. RESULTS AND DISCUSSINS 

Unsteady MHD flow through a horizontal channel 
consisting of two immiscible fluids between to non-conducting 
plate is studied. Separate solutions for each fluid are obtained and 
these solutions are matched at the interface using suitable 
matching conditions. The partial differential equations governing 
the flow are transformed into ordinary differential equations and 
solved analytically considering cosine, sinusoidal and exponential 
function of at . The effect of governing parameter such as 
viscosity ratio, periodic frequency parameter cot , Hartmann 
number presented graphically for cosine, sinusoidal and 
exponential function of cot . The parameters are fixed except the 

varying one, cot = 45° , a = h = P l0 = Pn= 1, and M=2. 

The effect of viscosity ratio a on velocity is shown in 
Fig. 1 for cosine function of at . As the viscosity ratio increases 
the flow is suppress in both the regions, the magnitude of 
suppression is large in region-I compare to region-II. Similar 
effects are observed for sinusoidal and exponential function of 
at in Fig. 4 and Fig. 7 respectively. 



M 
P 
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NOMENCLATURE 

Hartmann number 

pressure 

time 

average velocity 



Greek letters 

p fluid density 

/j viscosity of fluid 

a ratio of viscosity 

a frequency parameter 

cot periodic frequency parameter 

Subscripts 

1,2 quantities for region-I and region-II respectively. 
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Fig.1 Velocity profiles for different values of viscosity ratio a 



The effect of frequency parameter COt on velocity 
profile is shown in Fig. 2 for cosine, Fig. 5 for sinusoidal and Fig. 8 
for exponential. As at increases the flow decreases in Fig.2 and 
Fig. 8. For at = 0° and 180° display the velocity profiles for 
steady case as seen in Fig. 5. 

Fig. 3, Fig. 6 and Fig. 9 are shown the effect of 
Hartmann numbers on velocity profile for cosine, sinusoidal and 
exponential function of at respectively. The magnetic field will 
falttern the velocity distribution profile and reduced the flow 
which is the usual Hartmann flow problem. 
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Fig.2 Velocity profiles for different values of periodic frequency parameter ©t 
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Fig. 7 Velocity profiles for different values of viscosity ratio a 




Fig.8 Velocity profiles for different values of periodic frequency parameter cot 
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Fig. 9 Velocity profiles for different values of Hartmann number M 
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